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Abstract 

We construct approximate kink solutions of supersymmetric open string field theory 
at lowest level when non-local operators in the tachyon effective action are fully 
taken into account. To this purpose we derive two duplication formula? for products 
of incomplete gamma functions, which determine the level of approximation of the 
solutions. The time kink is an instanton of the Euclidean theory with a well-defined 
tunneling probability. The spatial kink solution represents an unstable D9-brane 
decaying into a stable D8-brane. We calculate the ratio of the brane tensions, which 
almost exactly reproduces the expected value. 
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1 Introduction 



The prototype of instanton in local scalar theories is the classical Euclidean so- 
lution for a double-well potential, — V(<&) ~ $ 2 — $ 4 , representing a scalar field 
rolling between the local maxima of — V. The study of instantonic solutions 
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is an essential tool to understand the vacuum structure of the correspond- 
ing Lorentzian theory with the potential upside down. In particular, one can 
understand whether (i) the ground state is degenerate or not, (ii) tunnel ef- 
fects can take place between the two localized minima of the potential V, (iii) 
such a penetration through the barrier removes the degeneracy of the ground 
state and what is the splitting, and (iv) the system oscillates between the two 
minima and what is the frequency of the resonant state [1-3]. 

The same problem is neither trivial nor of mere academic interest as far as 
non-local theories are concerned. In fact, the simplest example of a non-local 
scalar with a static (i.e., zero-momentum) double-well potential is provided by 
the tachyon of open string field theory (OSFT, reviewed in [4-7]). In the case 
of nonchiral bilocal supersymmetric OSFT [8-10], the Euclidean (1/2, l)-level 
effective equation for a homogeneous tachyon is (up to a field redefinition) 



where t is Wick- rotated time and m 2 < 0, s > and a > are parameters fixed 
by the theory. In the appendix we briefly review how to obtain this equation 
from the fundamental OSFT action. Solving Eq. (1), even numerically, has 
been proved to be a formidable task and a simplified version of the dynamics 
was proposed in [10], where the potential is purely local: 



While for local systems a classical rolling solution between two maxima is a 
textbook exercise for the double-well potential, non-local dynamics can (and 
in fact does) spoil the neat correspondence between classical rolling and kink 
solutions. In fact, there exist Lorentzian solutions which start from a local 
maximum and apparently 'climb' the potential up and down indefinitely [11]. 
The same phenomenon was previously noted also for the bosonic cubic poten- 
tial [12-16]. 1 

For certain values of the parameters, a kink solution of Eq. (2) has long been 
established numerically [21] (also on a cosmological background [22,23]) and 
by an existence theorem [24] . 2 3 On the other hand, little is known about the 
solutions of Eq. (1). The first goal of this paper is to find an analytic, bounded, 



1 Analytic background-independent solutions for marginal deformations of 
Berkovits' string field theory were found in [17-20]. 

2 Kink solutions of the p-adic equation are discussed in [12,21,25-28]; existence 
(and uniqueness) theorems can be found in [27-29]. 

3 The existence of a kink was also argued in [30] on Minkowski and in [31,32] 
on a cosmological background. However, the action was truncated at second order 
in time derivatives and, due to the fundamental difference between non-local and 
finite-order theories, this type of evidence could not be considered as conclusive. 



(d 2 t -m 2 )e- sd ^ = a^e sd "e , 



(1) 



m 2 )e- sd ^ = a& . 



(2) 
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(very- well) approximated instanton solution of Eq. (2). In doing so we shall 
provide a method which is viable also for the original system Eq. (1). As a 
byproduct, we are able to verify quantitatively how well the approximation 
of [10] works for kink solutions. In addition, all the basic physical quantities 
associated with the instanton solutions in the local double-well potential can 
be reproduced also in this non-local context: energy splitting, penetration 
barrier, resonance frequency and so on. 

Moreover, the time-independent kink profile (i.e., where t is a spatial coor- 
dinate) is a Lorentzian soliton describing the decay of an unstable Dp-brane 
into a stable D(p— l)-brane. This is Sen's first and second conjectures, which 
we verify explicitly by calculating the ratio of the brane tensions. The result is 
very close (~ 99.8%) to the theoretical value T p _i/T p = \^2ir (in a' units; we 
will explain this 'descent relation' [6] in Section 3.2). This datum adds to the 
evidence that the diffusion equation is a natural consequence of the conformal 
symmetry of string theory, as will be discussed in [33]. 

In [34] we constructed an explicit, global, approximate solution of the ef- 
fective equation of motion (2). We made use of the diffusion equation ap- 
proach developed in [11,35,36] (see also [13,21,22,26,37-39] in relation with the 
same method and [40,41] for another recent treatment of non-local scalars). 
Non-local systems are reinterpreted as living at the fixed point of a higher- 
dimensional flow along the diffusion time r. The scalar $(t,r) obeys a heat 
equation which localizes all non-local operators, so that they act as shifts along 
r. Choosing a sign distribution as initial condition for the diffusing system, 
the diffusion equation evolves it to the error function 



which solves Eq. (2) with very good accuracy. 

The problem in this and other non-local systems is that, in general, powers of 
the candidate solution do not obey the diffusion equation and the potential 
does not translate naturally along r. This is true also for the error function. Its 
validity and the level of approximation were checked locally via an asymptotic 
expansion and a matching of the leading coefficients [34]. Here we propose an 
analytic method which allows one not only to solve globally the same problem 
in a more efficient way, but also to face the much more difficult problem of 
nested non-local operators, as in Eq. (1). 

As a first step, which is also of mathematical interest, we shall derive an ap- 
proximate formula expressing products of incomplete gamma functions 7(01, z 2 ) 
X7(/9, z 2 ) ... as an incomplete 7 with rescaled arguments. The error function 
is the special case a — 1/2, which is deformed to general incomplete 7's by 
non-local operators. 




(3) 
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While in previous works approximate solutions were seeked for particular prob- 
lems, here we operate the other way round and ask to which systems, through 
the approximate scaling formula, the incomplete gamma function is an ana- 
lytic global solution. The answer not only includes the instantonic solution 4 of 
Eq. (2), as expected from [34], but also the highly non-local Eq. (1), previously 
inaccessible by analytical methods, and two cosmological solutions for, respec- 
tively, a cubic and exponential potential. We shall report on the cosmological 
cases in a companion paper [42]. 

The paper is organized as follows. In Section 2 we derive two approximate 
scaling formulae for incomplete gamma functions and discuss their relevance for 
diffusing systems. The instanton of susy OSFT for the approximate potential 
and its associated tunneling probability are discussed in Section 3.1, while in 
Section 3.2 we compute the brane tension of the Lorentzian soliton. Section 
4 deals with the exact non-local potential. The main results are summarized 
and discussed in Section 5. We shall work in units of the Regge slope, a' = 1. 



2 Duplication formulae and diffusion equation 



We shall be interested in non-local systems described by the scalar field 

. = . (o ,,) = 2^), (4) 
r(a) 

where 7 is the lower incomplete gamma function (see Section 8.35 of [43]): 

l(a,z 2 ) = 2 [ Z drr 2a - 1 e- T2 . (5) 
Jo 

Here z is an evolution variable we shall keep positive for simplicity, z = is a 
branch point for <3> except when a is half-integer or integer, where $ admits 
the odd (or even) analytic continuation to the whole real axis. In such cases 
$ is analytic at the origin. More explicitly, for any non-negative integer N 



/ 1 \ 2 7V_1 Z 2j + 1 

$(iV + -,z)=erf(z)-e- 2 V — — , (6) 

V 2' J K J ^ r(j + 3/2)' 1 ; 

2 N z 2j 

*(N + l,z) = l-e- £ — , (7) 

3=0 



4 When considering homogeneous configurations, the spatial volume is factorized 
from the Euclidean action and the following kink- type solutions are 'quantum- 
mechanical' instantons. 



4 



where the first sum is omitted when N = 0. Some properties we will make use 
of are the asymptotic expansions 



~2a 
z^O Z 



$(a z) ~ - (8) 

1 ' ' r(a + l)' { ) 

$(a,z) 2 ~°° 1, (9) 
the shift property 

z 2a e -z 2 

$(a + l,z) = $(a, z) - , (10) 

1 (a + 1) 



and 



o 7 2o-l - 

d-Ma,z)= , (11) 

<9 2 2 $(a, z) = <9,$(a, z) (^-^ - 2z) . (12) 



2. 1 Fzrsi duplication formula 



The problem is to conveniently rewrite powers of $ in a form which manifestly 
obeys the diffusion equation. We begin with a product of two $'s with different 
first argument, which is a double integral: 

S(a, *)<&(/?, *) = r( J m \l d ^ r^-Vf ^e-^ 2 ^ 2 ) . (13) 

In polar coordinates (f?, 6) the integrand function factors into a radial and 
an angular part. This integral is not easier to perform, since the integration 
domain, a square of side z with one corner at the origin of the {t\,T2) plane, 
is not isotropic. A useful approximation, adopted also in Debye's model of 
phonons, consists in replacing the square with the quarter of a disk centered 
at the origin and with radius R = Rz, where 1 < R < \/2. The first quadrant 
is q e [0,R], 9 e [0,tt/2]. Then, 

r(a)r(p) Jo Jo 

1 Ma + {3,R 2 ), 



giving the first duplication formula of this paper: 

$(a,z)$(P,z) m§(a + P,Rz). (14) 
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The constant R can be determined by choosing R so that to preserve the 
isotropic integration area, z 2 = ttR 2 /4. Hence R = 2/ ' y/n ~ 1.128 for any 
a and (3 (half-integer if Eq. (14) is to be considered on the whole real axis). 
Actually, a more accurate estimate of R is obtained by matching the behaviour 
of Eq. (14) at the origin, yielding 



k = 



r(a + i)r(/3 + i) 



2(<*+/3) 

(15) 



One can check that Eq. (14) is satisfied with good accuracy when taking 
Eq. (15). For instance, one can estimate the quantity 



A max = sup A (z) = sup 



LHS - RHS 



LHS + RHS 



(16) 



where LHS and RHS are, respectively, the left- and right-hand side of Eq. (14). 
Since we shall make use of the first duplication formula with a = 1/2 and 
(5 — 1, let us take these values. For the rescaling coefficient Eq. (15), 

/3\V3 

K=(^) -1.145, (17) 

and A max 1.3%, while using the isotropic value R = 2/y^ the maximum 
error is above 2%. From now on we shall ignore the latter. 

At this point a numerical method to improve the result is ready at hand. In 
fact, one can find the value of R minimizing the global error: 



5 = mi J dzA 2 (z), (18) 



'o 



where z = O(10) is sufficient. We get 5 w 1.0% and A max « 0.97% at R 
1.138, which is very close to the theoretical value (see Fig. 1). 



num 



An alternative criterion is to minimize the L 2 norm of Eq. (14) and integrate 
only (LHS — RHS) 2 in Eq. (18). In general this yields a higher A max but a 
similar result as far as the global accuracy of the solution is concerned. 

Since non-locality has been all encoded in argument rescalings of the scalar 
field, we will have to deal only with algebraic equations. Therefore Eq. (16) 
or (18) or the L 2 norm are sufficient to check the accuracy of the equations. 
Nevertheless, one can show that Sobolev-type norms of the form, e.g., W n = 
J dz[LHS^ — RHS*-™- 1 ] 2 , where (n) denotes the nth derivative with respect 
to z, recover about the same values of k obtained with Eqs. (16) and (18). 
For instance, we have checked that the first Sobolev norms Wi, W2, ...for 
Eq. (14) yield k's between the isotropic value 1.128 and Eq. (17), with errors 
comparable with those above. This indicates that the estimate of the error is 
reasonably 'stable' with respect to derivatives. 
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Fig. 1. A of the first duplication formula, Eq. (14), with k = (3/2) 1 / 3 (theoret- 
ical value, thick line) and R = R nnm ~ 1.138 (numerical value, thin line). Here 
(a, (3) = (1/2,1). 

2.2 Second duplication formula 



Taking a = (3 and iterating Eq. (14), one obtains the scaling formula 

{$(a,z)} n pa $(na,Kz) , (19) 
where 

[T(na + l)} 1/(2na) 



k = K\n, a) 



[T(a + 1)] 



l/(2a) 



(20) 



The value a = 1/2 is of particular interest, because it corresponds to the kink 
of open string field theory: 



In this case, 



erf (z) . 



r.2 + i 



(21) 



1/n 



(22) 



We can check that the error propagating from Eq. (14) stays very low. When 
(n, a) = (3,1/2), the maximum error is A max »s 1.8%, which is rather good 
considering we have resorted only to analytic tools. The minimum S is 5 ~ 
1.3% at k w 1.23, very close to the theoretical value n = (6/7r) 1//3 w 1.24. The 
maximum error is A max pa 1.4%, slightly better than the one obtained by the 
analytic approximation (see Fig. 2). 

For the particular case n = 1 + 1/a, the sup-norm of A depends on a only 
very mildy. In fact, a change in a results only in a shift of the position of the 
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0.020 - 




Fig. 2. A with the theoretical (thick line) and numerical (thin line) value of k. Here 
(n,a) = (3,1/2). 

maximum to slightly higher z, while its height is almost unaffected. For this 
reason Fig. 2 is a reference also for other values of a when n — 1 + 1/a. 

An interesting possibility is the analytic continuation to real values of n. The 
scaling formula Eq. (19) is still a good approximation, with errors of the 
same order of magnitude of the ones previously discussed. Finally, it is worth 
mentioning that Eq. (19) seems to work well even for negative values of n, 
proving it is rather 'stable' against the variation of its parameters. We are 
confident that the range of applications of both duplication formulae will not 
be exhausted within this paper. 



2.3 Diffusion equation 



Let us consider a cosmological-type setting where all fields are homogeneous 
and live in a Z)-dimensional universe characterized by a Hubble parameter 
H(t), where t is proper time. Defining 

t = 2^/rz (23) 

and the d'Alembertian operator (on scalar fields; metric signature —,+,••• , +) 
□ = — <9 t 2 — (D — l)Hd t , the incomplete gamma function obeys the diffusion 
equation 

(□ + ^$^1=0, (24) 
provided 

1 -2« 1 

H = W^T)t- (25) 
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One gets Minkowski spacetime when a — 1/2. For any s, 



e sD $ a, = e~ sdr <5> a, — = = $ a, , . (26) 

The action of an exponential non-local operator on <3> is that of a simple shift 
along the extra direction r, but only if the d'Alembertian and derivatives of 
r commute, [□,<9 r ]$ = 0. This condition can be imposed only if the Hubble 
parameter does not depend on r. 

Comparing Eqs. (12) and (24), we infer an identity involving only first-order 
derivatives: (2rd r + td t )Q = 0. 



We conclude this section with another formula we shall often employ. On a flat 
background, $(1, z) = 1 — e~ z does not obey the diffusion equation. However, 
one can show that 

e -^/(l+4/3) 

^'■^'-TW (27) 

To get this result, one can use the series representation of the non-local opera- 
tor and recognize the generating function of even Hermite polynomials [44,45]. 
A simpler alternative is to Fourier transform $. 



3 Supersymmetric OSFT: Approximate potential 



We are ready to write a non-local Klein-Gordon equation for the scalar field 
given in Eq. (4). We start with the simplest case where non- locality is con- 
fined to the kinetic term. The normalization of $ is fixed to 1 without loss 
of generality; this will only affect the normalization a of the self-interacting 
potential. The mass m of the field is tuned by r. 

Recasting the r derivative in the diffusion equation (24) in terms of Eq. (11) 
and making use of Eq. (10), one has 

(□-^)*(a,s) = -^&(a + l,*). (28) 
Let us define 



$ = e sD $ = $(a,Kz) , (29) 

,2 



2 (XK 2 oi i 

m = , a = —m , f 30) 

r 

where 
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k — k\1-\ — , a 

\ a . 



) 



(32) 



Taking the scaling formula (19) for n = 1 + 1/a, Eq. (28) for r — > r / k 2 becomes 



For natural n and a integer or half integer (solutions analytic at the origin), 
only (n, a) = (2, 1) or (n, a) = (3, 1/2) are possible. 

While in [11,34-36] the local limit $ = $ corresponded to a local system 
with the same type of potential, here it gives rise to an altogether different 
theory. This is because we are in the context of the inverse problem, where 
the dynamics is defined through the field and, in turn, the local limit of $ is 
trivial. 

The local limit corresponds to s — > 0, i.e., k — > 1 and, from the scaling formula, 
n — > 1 (and then a large). The right-hand side of Eq. (4) tends to zero at large 
a, in agreement with Eq. (33) in the same limit, $(a, z) ~ $(a + l,z), which 
admits $ = as the only solution. Conversely, when m 2 , a and n are kept 
fixed Eq. (33) formally reduces to the properly said local system. 

The case (n,a) = (3, 1/2) corresponds to the effective Euclidean equation of 
the supersymmetric tachyon at lowest level in the field expansion and where 
non-locality is transferred completely onto the kinetic term. For a homoge- 
neous field and □ = — d 2 , Eq. (33) becomes 



Since r is positive, Eq. (34) reproduces Eq. (2) as promised, where t plays the 
role of Wick-rotated time. Therefore, this kink-type solution can be regarded 
as an instanton (see Fig. 3). 

Let us fix r by requiring the mass to be the one of the string tachyon in units 
of the Regge slope: m 2 = —1/2. Hence, r = k 2 . Then, s = k 2 — 1 0.539 for 
the theoretical value of k, while 



for the numerical estimate of k. This is in excellent agreement with the value 
of s dictated by conformal symmetry in string field theory, 




(33) 



(d 2 -m 2 )e- sd ^ = a^ 3 . 



(34) 



s w 0.515 



(35) 



s = s* = 21n(3 3/2 /4) w 0.523. 



(36) 



We will come back to this point in the conclusions. 
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Since non-locality affects the kinetic and potential energy as a relative rescal- 
ing of the coordinate, the usual asymptotic properties of a scalar field are 
automatically recovered. Namely, at t — ±00 the effective potential domi- 
nates over the kinetic term and the scalar pressure is negative (P + E ~ 0, 
effective 'cosmological constant'), while near the origin t ~ the pressure is 
equal to the energy (P ~ E, stiff matter). 



3.1 Tunneling 



In this subsection we shall reproduce all the basic features associated with the 
non-local instanton $(2). Its Euclidean action (per unit spatial volume) is 



-at 



$(□ + m 2 )$ + -$ 4 + - 

— — 



(37) 



where □ = — df and the constant cr/4 has been added to set the zero energy 
of the potential at the minima $ = ±1. Variation of this action with respect 
to $ yields the equation of motion (34). Pulling back a perturbed classical 
solution $ = $ c i + 8§ into the action, where 5$ is a quantum fluctuation, one 
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can expand Eq. (37) as 



S E = S [^] + ^jdt6^(t) 



() ' Sk 6$(t) + .... (38) 



5<S>(t)5<S>(t) 



Two particular cases are of interest. The first corresponds to the trivial solution 
$ c i = ±1 describing the (Lorentzian) particle localized in one of the minima 
of the potential. Then, So[±l] = whereas the quantum fluctuation operator 
reads 

i 5 ^) = e sD (□ + m 2 + 3a) = e* D (-d 2 + c 2 ) , (39) 

\ / *cl=±l 

where u 2 = m 2 + 3a. If the particle were localized in one of the minima, uj 
would be the proper frequency of the oscillator and u/2 the energy of the 
vacuum state (H — 1). 

The second case is when the classical solution is the instanton given in Eq. (3), 
^ci = ^inst- Such solution describes a (Euclidean) particle rolling between the 
two maxima of the inverted potential and corresponds, in Lorentzian signature, 
to the tunneling between the two vacua. Then, So^mst] = Ek + V can be 
computed and the integral showed to converge. The kinetic term is (integration 
on R) 



e k = ~\ J dt$i„st<9 2 <li nst 
1 



( , 7 , dz<S> inst (z)d 2 z <S> hlst (Kz) 

(40) 



yVr(l + k 2 ) 
while the potential is (m 2 = —a) 



4 
2 

a^/r 



V = - I It (1 -L''l- -h 

dz 
N 



mst^inst i ^ inst 
2 



4 



+ 2<Z> inst (z) [®inst(z) ~ $inst 



(«*)]} 



+ — - (Vl + k 2 - V2^) 



Here TV is the numerical factor 



/ 



dz 



z)-l 



1.1201, 



(41) 



(42) 



which may also be estimated analytically by using the scaling formula for 
a — 1/2 and n = 4 and n = 2: 



12 



Jdz [$l st (z) - l] 2 » J dz 



7(2, 

3tt 



= 7T 



2 9/4 



1.1603. 



-2 7 



(43) 



As expected, the error on the 'exact' result (which we shall use from now on) 
is of order 1%. Plugging in the values 



K 



1 + s* « 1.523, 



(44) 



as dictated by string theory, one gets So ~ 0.612. This fixes the first term 
in Eq. (38). Concerning the second term, we have the quantum fluctuation 
operator for the instanton 



'PS* 

5$2 



= e 



sU 



3>=*inst 



□ +m z + 3ae~ sD & 



inst 



(45) 



The last (non-local) term in Eq. (45) can be calculated with a further use of 
the second duplication formula, Eqs. (19) and (20) with n = 2 and a = 1/2, 
together with Eq. (27): 



\ irr J 



= 1 



/ ]T t f 1 P 7rr+4s 



y/nr + As 



(46) 



Up to a non-local factor, the quantum fluctuation operator in the instanton 
background, Eq. (45), simply becomes the Schrodinger operator with an at- 
tractive Gaussian potential: 



'5 2 S E 

£$2 



-dl + J 



\pKT + 4s 



(47) 



Standard semi-classical manipulations show that all the physical effects in- 
duced by instantons are encoded in the quantity 

n = Ke' So[ ^ ] , (48) 

where 

V 2tt 

and, as usual, the prime over determinant means omission of the zero mode. 
The physical meaning of fl is that of the resonant frequency between the 



det (5 2 S E /5<S> 2 ) ±1 



— j-/^ 



(49) 
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two minima of the potential due to tunneling (and therefore the tunneling 
probability per unit of time). Moreover, twice its value is the energy splitting 
of the even and odd superpositions of the local ground states. 

The final step is then the evaluation of the ratio of determinants. To this 
purpose, the hypothesis of the Gel'fand-Yaglom theorem [46] are satisfied 
and we can trade in the calculation of functional determinants for a one- 
dimensional Cauchy problem. 5 

Since the potential induced by the instanton admits normalizable bound states, 
we can apply the Gel'fand-Yaglom theorem on the compact interval [-T, T] 
setting eventually T — > oo. For this part of the calculation, the non-local 
problem is very similar to the local one, where the Poschl-Teller potential is 
replaced by a Gaussian one. The ratio of the determinants is given by the ratio 
of the solutions of the corresponding homogeneous differential equations with 
Cauchy initial conditions, namely 



det' 



-d 2 t 



det [-d 2 t + u 2 } 



lim 

T->oo 



(50) 



where 



d 2 t 



_ t 

y/nr + 4s 



VwW = o, Vw(o) = o, Vw(o) = i. 



(51) 



(d 2 -co 2 )^ ±1 (t) = 0, V±i(0) = 0, ^±i(0) = l 



(52) 

The solution of Eq. (52) is immediate, ip±i(t) = smh(ut)/u, whereas the 
solution of Eq. (51) can be found numerically. To this purpose let us fix the 
parameters to the values suggested by string theory, Eq. (44), and remember 
that m 2 = —a, so that u 2 = 1o = 1. 6 Then, 



D, ^0.34. 



(53) 



This means that the barrier can be easily penetrated by the non-local system, 
and that the even superposition of the wave functions localized in the two 
minima has an energy eigenvalue much lower than that of a particle localized 
in one of the two minima: If -?/>l and i/j-r are, respectively, the ground states 



5 See [47] for a generalization to higher dimensions and [48] for an introduction. 

6 This is consistent with the WKB approximation used to define the transmission 
coefficient, requiring uj 3 3> cr/4; in our case 1 > 1/8, which is reasonable. 
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localized in the left and right minimum, with energy eigenvalue -El,r = ^/2, 
then the true ground state of the system is ip = (1/-\/2)('0l+'0r)> with energy 
eigenvalue E = u/2 — Q ^ 0.16c<j. 

With the same choice of couplings (— m 2 = a = 1/2), the instanton of the 
local system (Eq. (37) with r = and $ generic) is 

$iocai(0 = tanh(t/2) , (54) 

and one obtains 

^local ~ 0.65 , fil oca l « 0.33 . (55) 

Thus, concerning the tunneling frequency, the difference between the local and 
non-local cases is about 3%. Taking into account the level of approximation of 
the scaling formula and its propagation (that does not exceed 2% in the model 
under investigation) one can conclude that non-locality does affect tunneling 
events, although in a very mild way. These results are also compatible with the 
difference (in L 2 norm) between the local and non-local instanton solutions, 
tanh(t/2) and & inst (z): although small, it is well above 2%. Consequently, there 
is a quantitative departure from local physics. Conversely, it is remarkable that 
all the effect of non- locality is just a very mild modification of the physical 
parameters, especially in view of the drastic changes that non-locality provides, 
instead, in the Minkowskian classical solutions of the same problem. 



3.2 Brane tension 



In order to complete the proof that Eq. (3) is a solution of OSFT, we must 
verify that it describes a BPS brane configuration. We refer to [6] for a review 
on the subject. 

According to Sen's first and second 'conjectures' [49-54], 7 the tachyon vac- 
uum corresponds to an unstable (non-BPS) Dp-brane (whose volume V p+ i we 
normalize to 1), which decays into a stable configuration. Therefore, at its local 
maximum the effective tachyon potential equals the tension of the non-BPS 
Dp-brane, which is [49] 

T ' = *k ' (56 » 

where g a is the open string coupling. When p = 9, the brane coincides with 
the target spacetime of Type I/IIA theory. On the other hand, a tachyonic 



7 As a point of fact these are no longer conjectures, since they have been proven or 
strongly supported numerically in the level truncation scheme [55-59] and analyti- 
cally [60-65]. See [7] for an updated discussion on the subject and references. 
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kink solution $(x) interpolating between the two minima is interpreted as the 
BPS D(j> — l)-brane, localized at x — 0, into which the unstable Dp-brane has 
decayed [51-53]. The tension of this brane is 

Vi = V2nT p . (57) 

The prefactor takes into account reduction of dimensionality of the brane 
(Tp-i = 2irT p ) and the fact that the tension of an unstable Dp-brane is y/2 
times the tension of a BPS Dp-brane. 

If t — x is regarded as a space coordinate in the (p + l)-dimensional target 
spacetime, our solution (which is now a Lorentzian soliton) is a candidate 
realizing the above picture. To show this, we must first revert to the original 
effective action of string field theory (with approximated local potential) and 
then fix the normalization of the solution. The former is 



$ ( Ql + I) $ _ ^$4 _ 2A 
2 / 2 



(5S 



1 

where A is a constant and 

$ = e _sS *$ = $ erf (^j . (59) 

Here we used r — k 2 , r — s — 1. In the previous subsections we chose the 
normalization $ = 1. Otherwise, the quartic term is rescaled as 1/2 = a — > 
ct/$q = 1/(2$q). Now we can fix this residual free parameter so that V(0) — 
V(±%) = V(0) = T p , yielding 

<S> 2 = 8A = 8g 2 o T p . (60) 

The truncation level of the action affects the value the non-BPS brane tension 
T p and possibly the ratio 



(61) 

A non-trivial check to perform on the thus-established solution is to see 
whether Eq. (61) is satisfied. Explicitly, and using the second duplication 
formula, 



2tt 

^4.558, (62) 

where k = k(3, 1/2). This value reproduces the expected tension ratio y/2% w 
4.443 within 2.6%. Taking instead the string value Eq. (44), one gets —S*/T P 



p ~ 
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4.584, 103.2% the theoretical ratio. Evaluating the action exactly, i.e., without 
using the duplication formula?, one obtains 



SJ$ $1 r- 

■ 1 ' J « 4.435 = 0.998 x (V2tt) 



T -v.-, (63) 

for « = «(3, 1/2), and -S*/T p « 4.461 « 1.004 x (v^vr) for Eq. (44). 

Considering that erf(x/2) was regarded as the approximate solution of the 
lowest-level approximate effective action, the agreement is impressive. 



4 Supersymmetric OSFT: Non-local potential 



We now turn to the real effective (1/2, l)-level-truncated potential of OSFT 
which contains 'nested' non-local operators. In particular, we would like to 
reproduce Eq. (1), 

for some constants m 2 , s and a. Since the potential itself is dressed with a 
non-local operator, for purely technical reasons this system has been typically 
traded for its simplified form Eq. (33) in the literature. However, we are in 
a position to bypass this approximation and deal with the original case in a 
rather simple fashion. 

We can appreciate the advantage of using the duplication formula? when deal- 
ing with equations like (1). By construction, all the $ functions satisfy the 
heat equation and the action of the non-local operators is a translation in the 
evolution parameter r, Eq. (26). However, products of the scalar field generally 
do not satisfy the heat equation and therefore Eq. (26) cannot be used. The 
role of the duplication formulae is then clear: They select the best approxima- 
tion for non-local products of $ functions still satisfying the heat equation. 
Consequently, we shall reduce Eq. (1) to an algebraic equation, as was done 
also in [35] for cosmological models. 

The analysis of the previous section taught us that non-local effects are rather 
mild and the dynamics does not change much with respect to the local sys- 
tem. Therefore we do not expect the non-local operator on the left-hand side 
of Eq. (1) to modify qualitatively the behaviour of the solution with local 
potential. 

Indeed, an approximate solution of Eq. (1) is again the error function $(1/2, z). 
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Eq. (28) for a = 1/2 can be written as 



(64) 



On the other hand, according to Eq. (26) the non-local operator can be recast 
as 



2' y/T+W ' 



(65) 



Let m 2 = —1/2 and 
s = 5— r, 

AC 



(66) 



where and r are constants later to be determined. The right-hand side of 
Eq. (1) is a times 



'I' -) v-"'<\>- (J, : ~ 



' ! - «I> I -, - J <> <I>( 1. /,: j 



i_ _L e -(£*)' 



(67) 



while the left-hand side is l/(2r) times 



(^ + r)e-*«*(i. 



where 



(64) 2 
— rb 



(68) 



K = 



1 



K 2 — kI 



r — s 



(69) 



Incidentally, notice that a plot of <£> 3 against the right-hand side of Eq. (67) 
clearly shows that the maximum deviation from the local and the non-local 
case is at and near the origin, a fact which had been already realized in the 
comparison between local and non-local kinks with local potential [12,30]. 

Combining Eqs. (67) and (68) one obtains the algebraic equation 



/3 _ \ 2ra (\ 
* 2' fo ) + ^* 2' 21 



2m /l 
^* 2' 2 



£-1) *(!,«). (70) 
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At large z the correct asymptotic behaviour is obtained when a = 1/2. The 
simplest possible case is given by 

k* = k = k = -^= . (71) 

Then, one can check that 

A max « 0.6% , 5 w 0.7% , for r « 1.418 , (72) 

which is well below the error of Eq. (19). If one optimizes also for k, or employs 
the first duplication formula (14) in the left-hand side of Eq. (70) with R given 
by Eq. (17), A max is about the same. Amusingly, if instead one takes the 
isotropic value R = k, the maximum error drops below 0.04%! At any rate, 
the global accuracy of the solution is effectively determined only by the one 
of the duplication formula. 

The exponent s associated with Eq. (72) is s pa 0.304, which is about 40% 
off the string value Eq. (36). This should not discourage the reader for two 
reasons. First, it is not difficult to see that Eq. (3) is still a good solution 
of the system when Eq. (36) is enforced. In that = \Jl + K 2 s*/r 

and the integral in Eq. (18) is minimized with respect to r when r pa 1.620, 
for which A max pa 4.7%. One can refine the approximation by letting k free 
and fix the pair k, r so that the error on the second duplication formula for 
(n, a) = (2, 1/2) and the one on Eq. (70) are the same. The result is k pa 1.172, 
r pa 1.616, A max pa 3.8%. To improve this solution one could deform it locally, 
for instance as done in [26], with the requirement that the perturbed solution 
still obeys the diffusion equation. We shall not pursue this possibility here. 

Secondly, the brane tension ratio is about the same as, if not better than, the 
value obtained in the previous section. Take the Lorentzian string action with 
parameters given in Eqs. (71) and (72), 



S* = ± J d^x [$ + (e^$f - 2A 

where now r — s ^ 1. Using the second duplication formula one obtains 



(73) 



%> \J^ij' — s) 



+ 7T 



2 v / f^s-2 v /f+ v " ' 



^2(irr + As) 



w 4.517= 1.017 x (V2tt). (74) 

We conjecture that the exact result is extremely close to the theoretical value, 
although one should employ a careful numerical procedure (for instance the 
iterative method of [37,38]) to evaluate the fully non-local action. 
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5 Conclusions 



We can summarize the main achievements as follows: 

• We have derived approximate duplication formula? which express products 
and powers of incomplete gamma functions as a single 7 with rescaled ar- 
guments. Their level of approximation (about 1%) is under control and 
sufficient to easily construct instantonic tachyon solutions of supersymmet- 
ric string field theory, where both the effective kinetic term and potential 
are highly non-local. 

• The solution of the Euclidean OSFT problem at lowest-order truncation 
level with local potential is the error function erf (z). Despite the non- locality 
of the tachyon kinetic term has been fully taken into account, with respect to 
the 'local' instanton there is little quantitative and no qualitative difference 
in the physics of tunneling events. 

• The solution of the exact problem is still the error function, with the same 
global error. With respect to the solution for the local potential, the only 
difference is in the value of the parameters in the equation of motion. 

• To the best of our knowledge, this is the first time a solution of the non-local 
effective equation (1) is found either analytically or numerically. 

• Interpreted as a spatial kink, the solution in both versions of the theory is 
a Lorentzian soliton and, according to Sen's 'second conjecture', describes 
the decay of a non-BPS D9-brane into a codimension-1 BPS brane. This 
interpretation is fully confirmed by the brane tension ratio we obtained, 
which is in excellent agreement (99.8%) with the expected value. The role 
of non-locality was essential in the recovery of this result. 8 

Concerning the last point, the importance of lowest-level and non-local solu- 
tions has become apparent. Due to technical difficulties, off-shell calculations 
have been carried out only at lowest truncation level for both the bosonic and 
super symmetric string, so that the main effort has been focussed on the effec- 
tive spacetime action with non-local operators for the tachyon field only. This 
is sufficient to describe non-perturbative processes, as the tachyon is the only 
particle field with negative mass, and it is expected to encode the physics of 
unstable geometric configurations decaying into different others. Higher-level 
calculations as well as the recent analytic solutions in the full theory have con- 
firmed that the lowest-level tachyon dynamics captures most of the features 
of brane decays. Eq. (63) confirms this statement but adds to it the desirabil- 
ity of taking non- locality fully into account: Lowest-order calculations which 
are on-shell (i.e., local) or quasi- local (i.e., truncating non-local operators) are 
unable to achieve this degree of accuracy. 



8 For instance, taking the formal 'local' limit k — > 1 (s = 0) in Eq. (62) one obtains 
62% of the brane tension ratio. 
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We also notice that, for us, the tension of the single branes is unknown and 
is actually fixed by hand, Eq. (60). The ratio of brane tensions, however, is 
close to the expected value. In other words, non-locality can tell the difference 
between brane configurations with very good accuracy, but is insensitive of the 
properties of the individual branes. This is expected, as non-local operators 
govern the relative value of fields $(r, x) and $(r',x) in the equations of 
motion, but it does not affect the local structure of these equations, which is 
given by the level approximation. This is why the diffusion equation method 
plays an important role in the descriptions of the non-perturbative physics of 
very different in- and out-states. Its applications could be extended to other 
geometries (such as lumps of different codimension), but only after a spacetime 
non-local effective action has been constructed. This latter step is often very 
non-trivial and must be done with traditional techniques but, once this is 
accomplished, the rest would be under control. 9 

One of the most remarkable consequence of our procedure is the recovery 'out 
of the blue' of the parameters and the physics of OSFT, Eqs. (36) and (63), for 
the super symmetric local potential. We believe that the rather small difference 
between the actual values and those of OSFT is due to the approximation 
intrinsic of the duplication formula. This result is highly non-trivial. The mass 
and non-local exponent appear as separate inputs in the effective equation of 
the OSFT tachyon, although both are determined by conformal invariance. 
Obviously, a time rescaling can change their ratio, which is precisely the job 
done by the parameter r in our model. However, their product m 2 s = (k 2 — 
l)/2 is independent of r and is determined only by the value of k, which 
comes from the scaling formula for a particular type of potential. It attracts 
all our attention that the effective equation of the string tachyon with the 
same values of the coupling constants, as well as the brane descent relation in 
Sen's tachyon condensation, have been obtained starting from an apparently 
different framework. 

In other words, one can pose the question: What is the scalar field theory which 
describes tunneling between two inequivalent vacua and whose solutions obey 
the diffusion equation? The answer, surprisingly, is: The Euclidean tachyonic 
sector of open string field theory. The fact that string field theory may be 
viewed as a diffusing system was already pointed out in [11,34], where tachyon 
solutions of OSFT and boundary string field theory were mapped onto each 
other. We will formalize this picture in a forthcoming study and show that the 
diffusion equation implements conformal invariance at the level of the effective 
dynamics. 



Large-codimension lumps are expected to be described with a lesser degree of 
accuracy. Take the simplest case of a factorizable action (a typical example is the 
Gaussian case). The smallest error is obtained for codimension 1, because the prod- 
uct of more elements would propagate the error. 
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A Derivation of Eq. (1) from the OSFT action 

Contrary to the bosonic string, there are several proposals for superstring 
field theory, the first being Witten's [66-71]. The action was later modified by 
[8,72,73]. On a single non-BPS Dp-brane, 



where g Q is the open string coupling constant, / is the path integral over matter 
and ghost fields, Y-i is a double-step inverse picture-changing operator, a' is 
the Regge slope, Q is the BRST operator, * is a noncommutative product, 
and the string field ^± is a linear superposition of states (made of matter 
(super)fields X^, ^ and (super)ghosts 6, c, f3, 7) in the GSO(±) sectors [74], 
respectively. The coefficients of these states correspond to the particle fields 
fi of the string spectrum. 

The operator YL 2 can be either chiral and local [72,73] or nonchiral and bilocal 
[8] (see the literature and the review [4] for full details). These two theories 
predict the same tree-level on-shell amplitudes but different off-shell sectors. 
From now one we concentrate on the nonchiral version [8-10] and set a' — 1. 

Expanding the string field in an L eigenbasis, let hi be the conformal dimen- 
sion of the vertex operator associated with /j. One defines the level of fi to 
be M — hi + 1. By definition the tachyon has level 1/2. The truncation level 
(M, N) selects all quadratic and cubic terms of total level no greater than N 
made of fields with level < M [58,75]; in general, 2M < N < 3M [9]. 

In the picture and at level (1/2, 1), which is the lowest for the susy tachyon 
effective action, all particle fields in ty± are neglected except the tachyonic 
one, labeled $(#) and depending on the center-of-mass coordinate x of the 
string, and an auxiliary level —1 field u(x). The Fock-space expansion of the 
string field is truncated so that ^+(z) = j d p+1 k u(k)c(z)e lk ' x ^ and ^-(z) = 




(A.l) 
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J d p+1 k${k)-f{z)e ik - x( - z l The spacetime action reads [9,10] 




(A.2) 



where s* is given in Eq. (36), □ = d^d^ and <£> = e s * D//2 $. Combining the 
equations of motion for u and $, 



one obtains an equation of motion for $ alone [10]. It is the presence of the 
field u in the cubic action Eq. (A.2) which generates a quartic term in the 
tachyon effective action. 

Specializing to a one-dimensional spatial configuration (or a Wick-rotated ho- 
mogeneous background), Eq. (73) (respectively, Eq. (1)) is recovered for spe- 
cific values of m 2 , s and a. 
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